Fractionalization noise in edge channels of integer quantum Hall states 
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A theoretical calculation is presented of current noise which is due charge fractionalization, in 
two interacting edge channels in the integer quantum Hall state at filling factor v — 2. Because of 
the capacitive coupling between the channels, a tunneling event, in which an electron is transferred 
from a biased source lead to one of the two channels, generates propagating plasma mode excitations 
which carry fractional charges on the other edge channel. When these excitations impinge on a 
quantum point contact, they induce low-frequency current fluctuations with no net average current. 
A perturbative treatment in the weak tunneling regime yields analytical integral expressions for the 
noise as a function of the bias on the source. Asymptotic expressions of the noise in the limits of 
high and low bias are found. 

PACS numbers: 71.10.Pm,73.43.-f 



The fractionalization of the unit electron charge is an 
emergent phenomenon which occurs in a variety of low 
dimensional interacting electron systems [H-Q. The most 
known of these is the fractional quantum Hall effect [2J, 
in which the low-energy edge excitations carry a frac- 
tion of the unit charge that can be measured by shot- 
noise measurements [5|-|7[. A different kind of such frac- 
tionalization, which is the focus of this Letter, may 
also occur in the integer quantum Hall effect (IQHE) 
regime. Integer quantum Hall states that have filling 
factors (FFs) larger than 1 support several copropagat- 
ing edge channels [8(. If the electrons which flow in these 
edge channels are strongly coupled by Coulomb interac- 
tion, the edge excitations no longer have the usual Fermi- 
liquid-like behavior. Rather, the edge channels are de- 
scribed by the chiral Luttinger liquid theory, which pre- 
dicts that a single electron excitation in one edge channel 
separates into several copropagating plasma modes with 
different velocities. Each of these modes carry fractions 
of the unit electron charge in each of the channels, de- 
pending on the interaction strength. 

These fractional charge excitations in the IQHE have 
not been directly observed yet. However they may have 
had a crucial influence on recent experimental results; 
recently observed energy equilibration and energy loss 
in the electron transport at FF v — 2 9] suggests an 
energy transfer between the two channels without tun- 
neling. Controlled dephasing experiments of electronic 
interferometers [U 12 1 revealed a strong interchannel in- 
teraction at FF v = 2. In addition, charge fractionaliza- 
tion at FF v = 2 was raised as one of the explanations fl3j 
to the observed nontrivial behavior of the visibility of the 
Mach-Zehnder interferometer (MZI) as a function of the 



Figure 1: (color on-line) Measurement scheme. S,D and G 
denote source drain and ground, respectively. The green and 
the red lines represent the inner and outer edge channels in 
FF v — 2, respectively. A tunneling of an electron in the inner 
channel through QPC1 evolves into propagating slow and fast 
modes with fractional positive and negative charges on both 
channels, as represented by the small bulges above and below 
the channels lines 



frequency current fluctuations - both are zero. One may 
try to detect the fractional charges using high frequ ency 
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source bias. [14 1 



How can one measure these fractional excitations di- 
rectly? The basic idea would be to inject an electron 
to one channel through a tunnel barrier, and observe 
the fractional excitations on the adjacent channel. Note, 
however, that the fractional excitations affect neither the 
average current at the adjacent channel nor the low- 



measurement as was proposed by Berg et. al. 
which may be within reach with current technology 

In this Letter a different approach is taken, by consid- 
ering a mesoscopic device in which the fractional plasma 
modes impinge on a quantum point contact (QPC) which 
is placed on their way. The system is sketched in Fig. [1] 
The 2DEG bulk (the gray area) is assumed to be in a 
quantum Hall state with FF v = 2, with two coprop- 
agating edge channels on each edge of the 2DEG. We 
assume that no tunneling events occur between the two 
channels. The source contact marked "S" at the top left 
is biased by voltage V s relative to other contacts, and so 

2 

it injects a net electron current I s = \V S to each of the 
two outgoing edge channels, according to Landauer for- 
mula [18L 19]. Those electrons propagate to the right, and 
impinge on QPCO, which is tuned to selectively transmit 
fully only the outer channel to a grounded contact, and 
completely reflect the biased inner channel toward QPC1. 
QPC1 is tuned to allow small tunneling probability T\ of 
electrons in the inner channel from region I to region II. 
An extra electron which tunnels to region II in the inner 
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edge channel thus propagates along the top edge toward 
QPC2, presumably in a form of fractional plasma modes. 
QPC2, in turn, is tuned to allow small electron tunnel- 
ing probability T 2 in the outer channel from region II 
to region III. A tunneling event through QPC2 creates 
an extra charge on the outer edge channel in region III, 
which propagates to the contact at the top right of Fig. 1 
and adds temporarily to the current which is picked up 
by an amplifier. 

In the above setup, the fractionalization of the extra 
electrons in region II leads to low-frequency current noise 
in region III. Below, I first present a qualitative argument 
for the existence of this noise. The noise is then calcu- 
lated as a function the source bias V s at zero temperature 
in the weak tunneling limit of the QPCs, using the chiral 
Luttinger liquid theory. It is found that the noise behav- 
ior has a crossover; in the limit of large source bias the 
noise scales as V} +2ri where < n < 0.5 depends on the 
coupling between the channels (see definition below). In 
the low bias limit the noise vanishes as fast as V^, 3 . 

The appearance of the noise can be explained intu- 
itively as follows. Suppose that an electron tunnels to 
the inner channel in region II through QPC1. Because 
of the interaction between the two channels, the extra 
charge breaks up into two density modes, both propagat- 
ing to the right; a fast, "charge-like" mode, with negative 
fractional charge excitations on both edge channels, and 
a slow, "dipole-like" mode, with extra negative charge 
on the inner channel and extra positive charge (holes) 
on the outer channel (see Fig. [IJ . When each mode ar- 
rives to be near QPC2, it allows temporarily a certain 
tunneling event from region II to region III, in the outer 
edge channel. The charge-like mode induces a tunneling 
probability of an electron above the Fermi level, while the 
dipole-like mode induces a tunneling probability of a hole 
below the Fermi level. The two density modes have dif- 
ferent velocities, and so they arrive to QPC2 at different 
times, which are well resolved if the uncertainty in the 
energy of the tunneling electron is large enough (i.e. for 
high enough source bias). In this case the two corre- 
sponding tunneling events are separated in time and are 
statistically uncorrelated. As a result, the net extra cur- 
rent in region III is expected to be zero on average, with 
an equal average number of electrons and holes tunnel- 
ing through QPC2. However as the tunneling events of 
electrons and holes are random, the current will fluctuate 
around the zero average, and these fluctuations will have 
a low-frequency component, similar to the usual Schottky 
noise from a tunnel barrier. 

For a quantitative prediction for the noise, let us model 
the system by a low-energy effective theory in the lowest 
Landau level, using the Hamiltonian 

H = ^ Hp + Hqpci + Hqpc2- 

R 

The index R goes over the regions, R G {I, II, III}. Hp 



describes the evolution in the two edge channels within 
region R, corresponding to electrons with spin up and 
spin down relative to the magnetic field direction, with 
Coulomb interaction between the channels. In an appro- 
priate choice of gauge, one has (Ti = 1) 

dx Kn^in.R^^n.R + w °ut^ out ,R ^VWt,R 

u f°° f°° 
+— / dx : Pi n ,R/W,R : +5p,ieV s / dx : p in ,R(l) 

^ J — oo J-oo 

Here tpi n ( out), k( x it) is the electron annihilation oper- 
ator in region R in the inner (outer) channel (the 
Heisenberg picture is used throughout the Letter) and 
An(out).R(M) = ^i a (out),R^in(out),R are the Id electron 
number densities at the channels. ':' denotes normal 
ordering, -y; n and v on t are the bare velocities of the two 
channels and u is the interaction strength. The grounded 
contacts are modeled effectively by setting the integral 
boundaries to ±oo. The last term in Eq. ([1]) models the 
bias eVs of the inner channel in region I after QPCO. 

Hqpci describes the tunneling at QPC1, at x = 0, 
between the inner channels in regions I and II. Hqpc2 
describes the tunneling at QPC2, at x = L, between the 
outer channels of regions II and III. They are given by 



H 



QPCl 



Ci„V7\V>/ n „(0)Vvi(0)+H.c. 



(2) 



H, 



QPC2 = Wout 1 



^LiiiWVwjiW+H.c. (3) 



Here T\ and T 2 are the electron transmission probabilities 
of QPCl and QPC2, respectively. v- m and u out are the 
renormalized tunneling density of states in the inner and 
outer channels, which are assumed here for simplicity to 
be equal for all three regions. They cancel out in the 
calculation below and do not appear in the final formula 
for the noise. 

The measured noise in region III is given by (20j 



Sf-,o - 



dt (<S> eVs | {Iout,IIl(t), /out JIl(O)} |$ e y s 



(4) 

The state |$ e v 3 ) is the ground state of the unperturbed 
Hamiltonian Hp (i.e. with no tunneling events be- 
tween the various regions), where the inner channel of 
region I is biased by eV s and all other edge channels arc 
grounded. The current operators in Eq. (j4}, / ut,in = 
e l± dxp out ju, can be written using the tunneling op- 
erators, 



t 

out, III 



(■MW'out.nCM) - h.c. 



(5) 

The noise is calculated by expanding the time evolu- 
tion of the current operators to first order in each of 
the tunneling probabilities T\ and T 2 using Keldysh 
formalism [2 1[ |22j. 
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The evolution in region II is solved analytically. Note 
that given the Hamiltonian Hr = u in Eq. ([T]), The equa- 
tion of motion for the density is given by 



d_ ( Pin,II 
dt \ Pout, ii 



dx 



Pin, II 
Pout, II 



= 0, 



(0) 



where U= 



v m u 
u v out 

density modes of the dynamics are the eigenstates of 
U, which can be written in an orthonormal basis as 
cos 9 \ ( sin 9 

sin 9 J \ — cos ( 

velocities of the density modes are the eigenvalues of U, 



is the velocity matrix. The 



where < 9 < ir/2. The 



Vl,2 



l+^out 



± 



1 2 

-(Win-fout) +U 2 , (7) 



where the +(— ) sign refers to the fast(slow) mode. 

The noise in Eq. ((4]) is now calculated in two steps. 
First, we expand the time evolution of the current oper- 
ators to second order in ^/T^. Using Fourier transform 
to express the result as integral over energies, one finds 



Sf^o = 2g T 2 dujn L (uj)+2g T 2 du (1 - 71l(u))) , 

JO J -oo 

(8) 

2 

where go = tt is the unit conductance, and the function 
n^(w) is given by the correlator 



FT ' 



n L (oj) = v out ($ eVs | ipl ut U (L, 0)Vw,n(£, t) \$ eVs 

(9) 

where "FT" denotes Fourier transform, til(uj) can be 
interpreted as the effective mean occupation of electron 
states at energy uj in the outer channel in region II at 
x = L, near QPC2. Without inter-channel interaction, 
in the case u — 0, til(uj) would be a Fermi distribution 
at zero temperature, which is a step function tll(uj) = 
(— a;). However, when u ^ 0, the noisy inner channel 
excites the electrons in the outer channel and changes its 
occupation distribution function. One can write 



n L (uj) 



du) n 



2tt 



-<3 (-UJ + uw) B (wout) ■ 



(10) 



The function B (w out ) is related to the probability for 
an electron-hole excitation with energy u; ou t in the outer 
channel in region II. In the absence of net current in 
the outer channel, the electrons and the holes contribute 
equally to the noise. We can therefore consider only the 
first term in Eq. <j8j) twice, which leads to 



duj n 



2- 



B(uj 



^out-^V^out 



(11) 



The second step is to calculate the function B(ui ou t) up 
to second order in the transmission amplitude of QPC1, 



y/Ti, using the chiral Luttinger liquid theory (see sup- 
plementary material for details). One finds 



-B(aw) = TiT 



r eVs dWin / 
Jo 2vr \ 6 



eV s — U}[ n J C (cD ou t, Wi n ) 



(12) 

Here r = — — — is the relative delay of the arrival 
of the two modes from QPC1 to QPC2. We define the 
second mode to be the dipole-like mode, such that r > 0. 
In Eq. (|T2|) we also have eV s = eV a r, and w ut = u ou tT. 

The function C (cj out , o>i n ) weights the contribution of 
processes with energy loss Cj- m /T in the inner channel to 
the probability of electron-hole excitations with energy 
Wout/T in the outer channel. Note, however, that this 
function can have negative values. It does not depend on 
the bias, and is a property of the free evolution of the 
two modes from x — to x — L. It is found to be a sum 
of three terms, C (w out , w in ) = + + with 



C (2) 
C (3) 

C (4) = 



(2ir)\W(Q out )\ 2 5(Q in -ij out ) (13) 
23? {e^-^+^W (w, n ) W (w fa - u> out ) W* (£ QUt )} 
^ la +o out ) f <^L e ^2u w w * r-. n _ 

J 27T 

xIF(-Wout+w)VF*(aj in -w + w ou t), (14) 



where W (uj) = w(t)pT, and 

The power rj is related to the density modes of Eq. ([6]), 

i] = cos 9 sin 9 = 4m/ \J (v\ — v 2 ) 2 + 4it 2 . Thus, the func- 
tion W (oj) is directly related to the fractionalization ef- 
fect. Fig. ([2]) shows the function W (uj) for three possible 
values of the power rj. Note that it vanishes at negative 
values, and satisfies VF(0 + ) = 2niri. Also note that the 
power-law tail of W (uj) at high energies is directly re- 
lated to the power-law divergence of w(t) at t = 0. The 



r(i+»?) 1 



, where 



asymptotic behavior is W (Cj 1 ) 
r is the gamma function. 

The fractionalization noise S7->0) calculated according 
to Eqs. (11)-(15) , is plotted in Fig. [3] as a function of 
eV s . The contribution of the noise only from is also 
plotted in Fig. [3] In the limit of small bias, eV s <C 1, the 
leading order of the noise comes from the contribution 
of at small values of a)i n , that is from the value of 
\W(0+) | 2 . Equations CCE]) and ([T2j) then lead to 



5' 



/->olv- a «i 



47T77" 



~ \ 3 



(16) 



In the opposite limit, of eV s ^> 1, The noise is a power- 
law function of the bias. One can see from Fig.|3]that still 
the main contribution to the noise comes from C {2 \ In 
this limit the noise is dominated by the tail of \ W (<D ou t) | 2 
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Figure 2: (color online) W(£j), absolute value and phase, for 
r\ = 0.5 (magenta), r\ — 0.3 (blue) and i) = 0.05 (black) 



cent experimental results 

HEUd- Ifthe y 

are indeed 

related to the same fractionalization effect, then the en- 
ergy scale is at the order w 10 fieV for a device with 
typical length of 10 [im. The crossover voltage should 
be different for devices with different typical lengths and 
different edge profiles. Finally, it should be mentioned 
that the effect of finite temperature and the effect of the 
disorder in the edge channels in region II on the behavior 
of the noise were not discussed in the model above and 
deserves future study. 

I acknowledge B. I. Halperin, G. Viola, Y. Oreg and 
E. Berg for very useful discussions. This work was sup- 
ported by NSF grant DMR-0906475. 




Figure 3: (color on-line) The fractionalization noise, 
Sj^o/(2goTiT2T~ 1 ) as a function of eV s = eV s r, for r] — 
0.5, 0.3, 0.05 (magenta, blue and black solid lines, respec- 
tively). The corresponding dashed lines show the contribu- 
tion of the term C"' 2 ' to the noise for the same values of r\. 



at high o)i n . Approximating C ~ C*- 2 -* and using again 
Eq. (JT3D in Eqs. ([11]) and ([12]), one finds for this limit 



S 



2g T 1 T 2 T~ 



(2r ? + l)[r(l + 7 ? )]- 



-2 n 



(17) 

In summary, the low-frequency noise due to the frac- 
tionalization effect was calculated in the integer quan- 
tum Hall effect at FF v = 2 at zero temperature. The 
noise is a result of the fractional density modes impinging 
on a QPC and inducing excess tunneling events of elec- 
trons and holes through the QPC to the outgoing leads. 
The fractionalization noise is found to vanish faster at 
the low source bias limit, where the time of arrival of 
the two density modes to QPC2 is no longer well re- 
solved. Thus, the bias in which the crossover occurs in 
the behavior of the noise corresponds to the difference 
in the arrival times of the modes from QPC1 to QPC2, 
eV s I crossover ~ Tvr . This crossover bias may be roughly 
estimated, based on energy scales which appeared in re- 
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Details of the perturbative calculation of the noise 
using the chiral Luttinger liquids theory 

The derivation from Eq. (4) to Eq. (8) in the main text 
is via perturbation expansion in the parameter \fT2~ at 
zero temperature, and is rather straightforward - see for 
example Ref. 16 for similar derivations. The derivation 



from Eq. (8) to Eqs. (12)-(15) in the main text is less 
obvious and involves some aspects of the chiral Luttinger 
liquid theory, as is elaborated here. 

Using Keldysh technique to expand tit, (lu) in Eq. (9) 
in the main text to second order in the transmission am- 
plitude y/Ti, one finds 



where 



n L (w) = (-w) + v m v out Ti 



G [e(h)] {o 



~2^ 



dWr 



2tt 



dtdt^g 1 ^ (t, t 1 ,t 2 )e- lwt - iWm(t2 - tl) 



(1) 



(2) 



and the two correlators g^ h "(t,ti,t2) are given by 

g^(tMM) = (*ev.|V'in > n(0 ) *2)^t,n(^*)V'2 u t > n(i.0)^ in (0 > ti) |<J> e yJ 
g^ h \tMM) = (* e v.|^ 1 n(0,ti)V'2L 1 n(^*)^2 u t > n(i,0)^ n (0,t 2 ) |* e v.) 

I 



(3) 



The superscript "0" on the ip operators denotes an evo- 
lution according to the free Hamiltonian H\\ in Eq. 1 in 
the main text. Note that all the other correlators ap- 
pearing the Keldysh perturbation calculation of nj, (cj), 
which have different ordering of the four i\) operators, 
vanish. This was checked for each such correlator by 
performing the integration over t in the FT in Eq. (9) in 
the main text using contour integration in the complex t 
plane, and rearranging the outcome. 

Only the integral term in the r.h.s. of Eq. (fTJ) con- 
tribute to the noise Sf^Q. The derivation of the integral 
of G^ (uj,uj in ) is given here in details. The integral of 
£?w {<jj,Lo in ) can be calculated by similar steps and gives 
no contribution at all. The physical reason for that is 
that at zero temperature we injects only electrons to re- 
gion II through QPC1, and not holes. 

Following Eqs. © and ©, the correlator g^(t,ti,t 2 ) 
is calculated using bosonization technique, by expressing 
the electron fields operators with the boson fields of the 
two chiral modes, 



"0° u(x t) ~ y e i[coa0cf>i(x— vit)+sin6cf>2(x-V2t)] 



%jp (x t) — \ _ e i [ sine <t>l(x-vit)-cos6<t>2(x-v 2 t)] /£\ 

where A is the momentum cutoff. It is further assumed 
here that due to the presence of the grounded ohmic con- 
tacts, any finite size effects which result from boundary 
conditions in the x axis can be neglected. The two boson 
fields thus satisfy the commutation relations 



[<j>i{x),^>j{x')] = TtiSijsiga(x - x') (5) 
and satisfy the diagonalized equations of motion, 

<Pi = v i~S—i * = 1)2, 

ox 

where v\ , v 2 are the velocities of the fast and slow modes 
which are given in Eq. (7) in the main text. The expres- 
sions in Eq. ^ is then inserted into Eq. ©, which is 
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then calculated by comparing its r.h.s to the l.h.s of the 
following identity: for any four operators of a boson field 
{4>(xj)},j — l,.-.,4 and a set of constant ai...,a4 with 
vanishing sum, a k — 0, one has 

(6) 

with K(xj -x k ) = ($ eVs | [cj>{xj) - 4>{x k )] <j){x k ) \$ e v,)- 
Eq. ^ is proved by repeatedly using the Baker- 
Hausdorff formula e lA e lB e i{A+B) e -|[A,s] 

four times, and then using the identity of 
boson fields (* eV ,| e*'£* |$ ey3 ) 
e -5<*ev s lE fc ^<M^)] 2 |*ev a >_ 

Equations © + (0]) leads to a product of two correla- 
tors similar to the one in the l.h.s. of Eq. ©, one for each 
boson mode, (f> = (f>\ and <fi = <p 2l with different parame- 
ters ai,..,a,4,: In the correlator for (f>i the parameters are 
a\ = — cos#, ci2 = — sin0, 03 = sin#, 04 = cos#. In the 
correlator for <fi 2 the parameters are a\ — — sin#, 02 = 
cos#, 03 = — cos 6, 04 = sin 8. For both fields, given 
Eq. ([S|), the correlator K satisfies ^e K ^ = jq^y- Ap- 
plying the identity in Eq. ([6]) results in a product over 
all possible pairing of bosons operators, which can be 
written as 

g^(t,h,h) =gZ ut (-t)g? n (t 2 -t 1 )c(t,t 1 ,t 2 ). (7) 



The r.h.s. of Eq. ([7]) is a product of three functions. 
The first two are the lesser and greater correlators, which 
result from the all the pairings of boson operators whose 
origin are two ip operators of the same edge channel (ei- 
ther the inner or the outer channel). They are given by 



g>(t) = -i($eyJC,n(0»*)C,n(°> )l $ eV s > = 
11 

2lTVi n t + iS 
9 < ut (t) = i($eV^ltM L >°Wout,Il(L,t)\$eV,) 



where v in = v^ os ' H 'v s 2 in '' ' 6 and v out = vf n * 6 e . One 
should note here that the above correlators g > and g < 
are also the ones that appear in the expression for the 
current through a biased QPC. Therefore Vi n and v ou t 
are identified with the renormalized tunneling densities 
of states at the QPCs which appear in Eqs. (2) and (3) 
in the main text. 

The function c(i, ti,t 2 ) in Eq. ([7]) results from the pair- 
ing of boson operators whose origin are two ip operators 
from different edge channels. There are four such pairing 
possibilities, hence one finds 



c(t,ti,t 2 ) 



1 + w ( J 



! + «;( 



t - U 



T 2 . 



1 + W( ) 



l + w*(- 



-h - t 2 , 



(9) 



where the function w (x) is given in Eq. (15) in the main text, and where t% = L/v\, t 2 — L/v 2 and r = t-i — T\. 
Note that for noninteracting case, w = and hence c(t,ti,t 2 ) = 1. 

Equations (fTJ), @ and ©-© leads to Equations (10)+(12) in the main text; when inserting equations ©-(0) into 
Eq. ([2]), performing the Fourier transform in Eq. © and inserting the result into the integral term in Eq. (fTJ), one 
gets equations (10) + (12) in the main text, with the function C ((!>„, a) out ) given by 



C (&in, &out) 



dtdtidt 2 c(t, ti, t 2 )e 



(10) 



Note that C (wi„,(Ii mi ) gets contribution from terms in 
the expansion of Eq. (J9j) that contains product of two, 
three and four functions w(x). These correspond to C^ 2 \ 
C<-V and in the main text. The contribution of terms 
containing only one w(x) vanishes. In this context, note 



that the calculation is greatly simplified by the fact that 
the Fourier transform of w(x) vanishes for negative fre- 
quencies, which causes the contributions to the noise of 
most of the terms in the expansion to vanish. 



